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Abstract 

^ | A scalar bimetric theory of gravity with a preferred reference frame 

i . is summarized. Dynamics is governed by an extension of Newton's sec- 

kjT), ond law. In the static case, geodesic motion is recovered together with 

Newton's attraction field. In the static spherical case, Schwarzschild's 
metric is found. Asymptotic schemes of post-Newtonian (PN) and 
post-Minkowskian (PM) approximation are built, each based on as- 
sociating a conceptual family of systems with the given system. At 
the 1PN approximation, there is no preferred-frame effect for photons, 
hence the standard predictions of GR for photons are got. At the OPM 
approximation, an isolated system loses energy by quadrupole radia- 
tion, without any monopole or dipole term. Inserting this loss into 
the Newtonian 2-body problem gives the Peters-Mathews coefficients 
of the theory. 

1 Introduction 

According to a common belief, a scalar theory should not be able to explain 
the effects of gravitation on light rays, which have been predicted on the basis 
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of general relativity (GR), and have received experimental confirmation. It is 
true that the most obvious covariant scalar theory does not predict any light 
bending pQ, for instance, but this particular example does not prove that 
no scalar theory can predict light bending. In the same vein, it is usually 
thought that a scalar theory should lead to monopole and dipole gravitational 
radiation, and that, therefore, the effects of gravitational radiation should be 
very important according to a such theory — contrary to observation. The 
latter point is believed to be even more true if a preferred-frame scalar the- 
ory is considered. For instance, Carlip 2j writes: "For a [scalar] non-Lorentz 
invariant theory, conservation laws place no restrictions on radiation, and 
monopole radiation should dominate, permitting aberration at all orders." 
However, the latter statement is in the form of a conjecture, which has to be 
checked for each specific theory. 

The aim of this paper is: i) to present to the readers of this Journal a 
summary of the scalar theory with a preferred reference frame jSJ HJ Gi] and 
of the proof that this theory predicts just the standard effects of GR on 
light rays; ii) for the first time, to investigate gravitational radiative energy 
loss in this scalar theory. In particular, it will be shown that this is indeed 
a loss, and that there is neither monopole nor dipole gravitational radiation 
for this theory. 

2 Summary of the scalar ether-theory 

The starting point of the theory is the heuristic interpretation of gravity as 
resulting from the pressure force exerted on matter at the scale of elementary 
particles by an imagined perfect fluid or "micro-ether". This interpretation 
leads to a definite set of equations: it is this set that really constitutes the 
theory, and that should be assessed from the predictions it leads to and the 
comparison of these predictions with experimental data. According to this 
heuristic view, gravity would be simply "Archimedes' thrust in the ether"; 
it turns out that Newtonian gravity had been interpreted in this way by 
Euler in 1746 jjj- The construction of the theory along this line is given 
in Ref. j^j, but the formulation of motion was later modified in order to 
obtain a conservation equation for the energy 1 Since then, the general 

1 The publication dates are confusing here: the first version of Ref. dates back to 
end 1992. 
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motion of a test particle is governed by an extension of Newton's second 
law (formerly [3], this was true only in the static case); this applies also to 
light-like particles (photons) jH]. The dynamical equations that apply to a 
continuous medium have been deduced from the dynamics of a test particle, 
by considering first a dust continuum These dynamical equations allow 
to define the gravitational extension of the Maxwell equations in the theory, 
which is consistent with photon dynamics as governed by the extension of 
Newton's second law (jH], part 3). 

2.1 Space-time metric and field equation 

The equations are given in the preferred reference frame E, which is the rest 
frame of the reference body M or "macro-ether". 2 The reference body M 
is assumed to be equipped with an Euclidean space metric g°. Observers 
bound to the frame E would measure the spatial distances defined with g°, 
and would also measure the "absolute time" T, if there were no gravitation. 
However, an "absolute" version of Einstein's equivalence principle implies 
that physical standards of space (resp. time) are contracted (resp. dilated) in 
a gravitational field jSj- The scalar field of the theory is the contraction factor 
(3 affecting the time interval dt^ measured by a clock fixed at point x bound 
to the frame E (due to the dilation of the clock period in the gravitational 
field), or the square / = (3 2 . The latter defines the 700 component of the 
physical, curved space-time metric 7, in any coordinates (x M ), with x° = cT 
(and c the velocity of light), that are bound to the frame E: 

dtJdT = y/f = /3<l. (1) 

The interpretation of gravity as Archimedes' thrust leads to a definite ex- 
pression of the acceleration vector g of gravitation, in terms of the gradient 
of the "ether pressure" p e , the gradient operator being defined with the phys- 
ical space metric g [3]. It turns out [I] that g may be reexpressed in a very 

2 A reference frame can be physically thought as a spatial network of observers equipped 
with measuring rods and clocks, allowing to define space-time coordinates (z p ) (jj, = 
0, ...,3) such that each observer has constant space coordinates z % (i — 1,2,3). Such 
coordinates are called "adapted to the frame" , see Cattaneo [§] • The 3-D reference body 
associated with a frame may be defined as the set of the world-lines of the observers, and 
may be described by the space coordinates (z 1 ). 
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simple form in terms of the Euclidean metric g° and the scalar /: 



g = -jVof, V / = grad g0 /, (grad^J? = g^fj. 



(2) 



(Matrix (g 0ij ) is the inverse of the component matrix (<7oy) of the Euclidean 
space metric g°.) Due to the "absolute" version of the equivalence principle, 
the field (3 defines also the dilation of an infinitesimal distance dl measured 
in the direction of vector g, as compared with the distance dl° evaluated 
using the Euclidean metric: dl = dl°/j3, while the infinitesimal distances 
in directions perpendicular to g are unaffected. This assumption may be 
reexpressed in saying that the physical space metric g in the frame E has the 
following relation to the Euclidean metric g°: 3 



(with g* = (gj), gj = g®jgi and g 2 = <?^g*g' J ). The space-time metric writes, 
in any coordinates (x^) bound to the frame E, with x° = cT: 



that is, 700 = /, lij = —gij,loi — 0. That the 7 j components are always zero 
in the ether frame E, is the natural translation of the assumption that this 

3 We note that Eq. J2J does not apply when g = 0. Indeed, the assumed contraction 
occurs in the direction of the gravity acceleration g, hence it is a priori undefined at those 
places where g = [unless the scalar potential / has the same value in that place as the 
one it has "at infinity" (i.e., / = 1): then, there is no contraction at all]. In the generic, 
physically realistic case, the vector g vanishes merely at isolated points, e.g. at one point 
between two massive bodies, and also at one point near the centre of each massive body. 
At an isolated point x such that g(x) = 0, the direction of contraction cannot be defined 
uniquely. However, Eq. applies in the neighborhood of x and shows that the metric g 
is discontinuous at this isolated point, but remains bounded in its neighborhood. Hence, 
geometrical elements such as the length of a line arc defined as with a continuous metric. 
Moreover, the examination of the case with spherical symmetry shows that no difficulty 
occurs in the dynamical equations either. A more difficult case occurs in the highly 
idealized situation of "collapse in free fall" : g vanishes in a bounded volume domain [3] . 
This case was easily solved in Ref. [3] for the relevant subcase of spherical symmetry, and 
the way of solution can be extended to the case without any symmetry. From a theoretical 
point of view, one should remind that, in this theory, the "true" geometry is in terms of 
the Euclidean metric g°, whereas g has merely the status of an "effective metric" — as is 
also the case in the "relativistic theory of gravitation" (RTG) |1U1 ITT] . 
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(3) 



ds 2 = r ) liv dx^dx u = f(dx ) 2 — dl 2 , dl 2 = gijdx l dx\ 



(4) 
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frame admits a global simultaneity defined with the absolute time T. 



In many problems, one may neglect the time variation of the "reference 
ether pressure" [3], pf(T), because the analysis of cosmo logical models [T2] 
shows that this variation takes place over long time scales, of the order of 
10 8 years. Then the field equation |3] may be rewritten 0] as: 

Af-Ulf] (,° = cT), (5) 



/V/7,0 c ' 2 

with A = diVgOgrad g o the usual Laplace operator, defined with the Euclidean 
metric g°, and where G is Newton's gravitation constant, and a = (T 00 )e 
is the mass-energy density in the ether frame. (T is the "mass tensor", i.e. 
the energy-momentum tensor in mass units.) This equation is merely space- 
covariant. It is immediate that, in any static situation, Eqs. (J2J) and 
associate just the Newtonian gravity acceleration with the given constant 
mass-energy density a = er(x). However, a = T 00 does not only involve the 
density of the pure matter energy (e.g., for "dust", the latter would be the 
rest-mass plus kinetic energy), for it is also increased by the gravitational 
field itself jl]. In the static spherically symmetric situation, the space-time 
metric (jlj) predicted from the field equation (jSJ and the boundary conditions 
(f(r) — + 1 as r — > oo and f(r) bounded as r — > 0) is quickly got in spherical 
coordinates r, 9 and (p for g°: 



2V(r)\ 2 2 1 2 2 7 2 r n 2 //)2 , -2/17,2 



da- 8 = 1 ^ c z c/T^ : - 7rT dr^ - d^, dQ z = d9 z + sin z 9 dip 1 

\ c z J l _ fLp. 

(6) 

where V(r) is the Newtonian potential associated with the density a = a(r): 
V(r) = r GM ( U > du, M(r)=47r f u 2 a{u)du. (7) 

Jr U JO 

(M = M(oo) = limr^oo M(r) is assumed finite.) Being a static solution, this 
solution assumes that the mass-energy density a is time-independent (and 
spherically symmetric) in the preferred frame, i.e., the massive body is at 
rest in the "ether" or in a stationary rotation with respect to it. If <r(r) = 
for r > R, then V(r) = GM/r for r > R, this is the Schwarzschild exterior 
solution. As to the Schwarzschild interior solution of GR, it is only known 
when, in addition, a(r) = Const, for < r < R, but Eq. © does not give 
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the Schwarzschild interior solution in that case. Note that, in contrast with 
the situation in GR, Eq. (JBJ) is valid as soon as M(oo) < oo and makes no 
difference between the "interior" and "exterior" cases. 

2.2 Dynamical equations 
2.2.1 Motion of a test particle 

This is governed by a natural extension of Newton's second law, the latter 
being defined in terms of the physical, curved metric. In the case of a constant 
gravitational field, it is relatively straightforward to extend Planck's special- 
relativistic form of Newton's second law (F = dP/dt with F the force and P 
the momentum) to the situation with gravitation, because the time- derivative 
of a vector (like P) must undoubtedly be defined as the "absolute derivative", 
the latter being induced from the covariant derivative with respect to the 
constant space metric g. In coordinates, the absolute derivative of a vector 
w depending on the parameter £ along a trajectory in a Riemannian space 
equipped with a fixed metric g is given by: 



with T l - k the Christoffel symbols of metric g. This derivative was used by 
Landau & Lifshitz ^3] to rewrite the spatial part of Einstein's geodesic equa- 
tion of motion as Newton's second law, for a constant gravitational field. In 
the case of a variable gravitational field, we still have, in the frame E (as also 
in any other possible reference frame), a well-defined space metric g — but 
now it depends on time, and we have a local time t x . We have found fTJ] 
that a unique definition may be given for the time-derivative Dw/D^ of a 
space vector w depending on the parameter £ along a trajectory in the given 
space M endowed with a metric g = g^ that varies with parameter £, under 
the following essential requirements: (i) it must be a space vector depending 
linearly on w, (ii) it must reduce to the absolute derivative D w / D£ if g 
does not actually depend on £, (iii) it must be expressed in terms of the 
space metric g and its derivatives, and (iv) it must obey Leibniz' derivation 
rule for a scalar product, i.e. 





d_ 



(S(w,z)) =g w, — ) + g 




Dw 



z 



) 



(9) 
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This unique definition is as follows |14j : 

Dw D w 1 _j 

WW + 2^ W ' ' = * - af (10) 

To define Newton's second law, the trajectory of the test particle is parametrized 
with the local time: £ = t x , synchronized along the trajectory as defined by 
Landau & Lifshitz [13J. Actually, definition (fTUj) is also valid (and compelling 
under the above requirements) for any reference frame in any theory with 
a curved space-time metric In the present case, the synchronized local 
time is simply obtained by integration of Eq. ([T} along the trajectory, be- 
cause the absolute time T is globally synchronized in the ether frame E (i.e., 
because the 7 i components are assumed zero in E). 

There is no ambiguity for the definition of the momentum, either. This 

is 

P = m(v)v, m(v) = m(v = 0)/y„ = m(0).(l - v 2 /c 2 y 1/2 , (11) 

the velocity v of the test particle (relative to the frame E) being measured 
with the local time i x , and its modulus v being defined with the space metric 
9 

* s K s jfW' ' s *^ = ^ (12) 

The gravitational force is assumed to be F g = m(v)g with g given by Eq. (J2J), 
and there may be also a non-gravitational force F . Thus, the extension of 
Newton's second law is obtained as 

DP 

F + m(t;)g = — . (13) 

For a photon, the velocity v is still defined by (fT2|) . but the mass content 
of the energy E = hu is substituted for the inertial mass, the frequency v 
(number n of oscillations per unit time) being measured with the local time: 

^ , dn , s 

E = hv, u= — , (14) 

E DP E 

F » + ? g = is;' p "? v - (15) 
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Defining the "pure energy" of a mass point (i.e., not accounting for the 
potential energy in the gravitational field) as 



E = m(v)c 2 , (16) 

we see that Eq. (j!5j) applies to both mass points and photons. Together with 
the assumed vector g (Eq. Q) and the assumed metric (Eqs. (j3J) and 
Eq. (JT5|l implies that free particles (F = 0) follow space-time geodesies in 
a static gravitational field (f } o = 0): see Ref. [H] for mass points and see 
Ref. [3] for photons. Also, Eq. (JT5|) with F = implies jl] a simple evolution 
equation for the energy E: 

<^ = *§ iP^f). (17) 

Using Eq. (|T7|). Newton's second law ([T3)l can be recast to give the 
"coordinate acceleration" for a free test particle (mass point or photon): 

^! = \ + 2/3 jU A ^ - r>% fc - V j ^« fc - - /(Vo/) 4 . (18) 
(3 \dT H ' J J 3k 2 y dT 2 J K 0J J K ' 



2.2.2 Continuum dynamics and energy conservation 

Equations (fTUj) (with £ = t x ) and (JT5j) (with F = 0) define Newton's sec- 
ond law in the preferred reference frame, for any free test particle. It is 
thereby defined also for dust, which is a continuum made of coherently mov- 
ing, non-interacting particles, each of which conserves its rest mass. Dust is 
characterized by its mass tensor 

T i*> = p *u»U u (19) 

(with p* the proper rest-mass density and f/ M = dx^/ds the 4- velocity). 
Independently of the assumed form for the space-time metric 7 (provided it 
satisfies 7oi = in the preferred frame), this extension of Newton's second 
law mathematically implies [Hj the following dynamical equation for dust: 



= K (20) 

where is defined by 

bo(T) = ±g jkfi Ti k , h(T) = -^g ik!0 T ok . (21) 
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(Indices are raised and lowered with metric 7, unless mentioned otherwise. 
Semicolon means covariant differentiation using the Christoffel connection 
associated with metric 7.) Equation (J2UJ) . with the definition (}2Tj) . is assumed 
to hold true for any material continuum: accounting for the mass-energy 
equivalence, this is the expression of the universality of gravitation. Equation 
(J2(Jj) is valid in any coordinates (x^) that are adapted to the frame E and 
such that a; = <f)(T) with T the absolute time. (In contrast, the validity 
of the scalar field equation (jSJ) demands that x° = cT.) For dust, Eq. ()20|) 
implies, conversely j^j, that the constitutive particles obey Newton's second 
law, and that the rest mass is conserved, i. e. 



{p*U L 



0. 



(22) 



For a general continuum, Eq. (J2*U|) implies an exact energy conservation, 
which substitutes [3] for the exact mass conservation (|22j) . With the help 
of the following relation, valid for the assumed space-time metric (Eqs. © 
and ©): 

-7 = 700-3 = f-(9°/f) =9° (7 = det(T^), g = det(^), g° = det(g^)), 

(23) 

one obtains from (J2U|l . in Cartesian coordinates (x l ) (i.e., = 5ij) and with 
x° = cT: 

Tl v = h T™= l -f fi o, (24) 



^i,v — 2 ° ~ 2 9ikfi T° k — ~ gjk,i T jk . 



(25) 



Using the scalar field equation fl^J), Eq. (J2l*j) may then be rewritten as [I] 



rpv , 



8ttG 



+ div (/ f 



.0 



(26) 



This is a true conservation equation for the quantity 

2" 



e = c 2 T ° + 



87TG 



2 , C 4 ( f,0 



(27) 



which may be called the total energy density, while e m = c 2 Tq and e g are 
the volume densities of the matter energy and the pure gravitational energy. 
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Integrating (J26|) in a fixed domain Q, and assuming that there is no flux of 
matter accross its boundary dQ, we get 



(V and S are the volume and surface measures associated with the Euclidean 
spatial metric g°] in (j26H28j) the scalar product and square also refer to g°). 

3 Effects of gravitation on light rays 

The prediction of gravitational effects on light rays in GR is done by assuming 
that the world lines of the light rays are null geodesies of the space-time 
metric, and that there is just one massive body, which is static and spherically 
symmetric (see e.g. refs. PUE3D-6J). This seems justified: in view of the 
(a posteriori) smallness of these effects, the deviation from a uniform motion 
becomes sensible only when the ray becomes very close to one massive body, 
so that the field of the other bodies may indeed be neglected; and the small 
deviation from sphericity is a second-order effect. As we have recalled above, 
in the static situation, the extension of Newton's second law assumed 
in the present theory implies f3j that photons follow space-time geodesies 
of metric 7, and these are indeed null geodesies, since the velocity v of the 
photons (as defined by Eq. (fT2*|) ) is the constant c. 4 Moreover, in the 
static spherically symmetric situation, the metric 7 © predicted by the 
present theory outside the body is just Schwarzschild's metric. One thus 
expects that this theory should predict the same effects on light rays as the 
standard effects predicted with the Schwarzschild metric in GR. However, 
in this preferred-frame theory, we cannot assume that the relevant body is 

4 In Ref. 0, photons are defined as light-like particles in the sense that v — c, in the 
same way as in GR. (Of course, v = c is equivalent to say that the tangent vector to the 
world line of the particle, U = (U^) with = dx^ /<i£, is an "isotropic" 4-vector, i.e. 
7(U,U) = 0.) The assumption v — c is compatible with Newton's second law (|15|) with 
Fo = . Indeed, the coordinate acceleration (|18(l can be rewritten to give 



(where g.v = t/(g, v)). By calculating d(v 2 )/dT, is easy to check that any solution of this 
with an initial data such that v{T = 0) = c satisfies v = c at any time. 



dE 




(28) 




(29) 
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at rest in the preferred reference frame E: we should account for a uniform 
translation of the body through this "ether" . 5 

3.1 Asymptotic post-Newtonian approximation 

To do this, we must develop an approximation scheme valid for weak gravita- 
tional fields and for low velocities of celestial bodies, thus a "post-Newtonian" 
(PN) approximation scheme. This should be done in following as closely as 
possible the rigorous procedures which are currently used in applied mathe- 
matics, because relativistic celestial mechanics aims at very high precisions. 
A PN approximation has been initiated along this line for GR by Futamase & 
Schutz [19J, and another one, which bears several similarities with the former, 
has been (independently) proposed and developed for the scalar ether-theory 
in Refs. jHJ 120] • We call this scheme the "asymptotic" PN scheme, because 
it is based on the usual method of asymptotic expansion, as applied to a 
boundary value problem for a system of partial differential equations. For 
light rays, we are in the simpler situation just described: a light-like test 
particle in the field of a spherical body that undergoes a uniform translation 
in the frame E, hence we need merely a summary of this scheme. 

3.1.1 Principle of the asymptotic PN approximation 

Since Eq. (J3J is a second-order hyperbolic equation, a gravitating system 
made of barotropic perfect fluids is defined (20] by the initial data for / and 
&rf (for the scalar gravitational field), and by the initial data for the pressure 
p and the absolute velocity u = dx/dT (for the matter fields). We want to 
describe gravitation in the Newtonian limit where: i) the space metric g (J3J) 

5 If the "absolute" velocity vb of the body is negligible with respect to c, we may consider 
that vb is indeed a constant, at least when calculating the trajectories of light rays. Now 
this theory predicts preferred-frame effects in celestial mechanics, and these have to be 
small, in view of the fact that Newtonian mechanics (which has no preferred-frame effect) 
is already very accurate. Hence it is necessary that vb be indeed small. The adjustment of 
the equations of motion for the mass centers |17| on observations will allow to determine a 
"best-fitting" value of the absolute velocity V of the mass center of the solar system, hence 
giving access to absolute velocities of all solar-system bodies. A preliminary adjustment 
(neglecting the effects of self-rotation of the celestial bodies, and assuming that the solar 
system is isolated) has led to the value |V| ~ 3 km/s ^H]. Hence one expects that the 
absolute velocities of solar-system bodies should have the same order of magnitude as their 
orbital velocities and, in particular, should indeed be negligible with respect to c. 
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becomes close to the Euclidean metric g° and the local time (JI} becomes 
close to the absolute time — both being true if / — > 1, and ii) the matter 
fields become close to corresponding Newtonian fields. To do that, we use 
an exact similarity transformation that exists in Newtonian gravity ^3 120] : 
If p^ l \ u^, = F^(p^) (density), and U-jp (Newtonian potential) obey 
the Euler-Newton equations, then, for any A > 0, the fields 

p w (x, T) = AV 1} (x, y/X T) , p< A > (x, T) = Ap (1) (x, T) , (30) 

(/f ) (x,T) = A[/; ) (x 1 VAT), uW(x,T) = v / AuW(x,VAT), (31) 

also obey these equations — provided the state equation for system Sa is 
pW^pW'j = \F^(X~ 2 p^'). This similarity transformation defines the weak- 
field limit in Newton's theory itself: as A — > 0, the potential and the density 
in the bodies decrease like A (while the bodies keep the same size), the (or- 
bital) velocities decrease like \/A, and accordingly the time scale increases 
like 1/y/X. Therefore, to define the Newtonian limit in the scalar theory, we 
just apply the similarity transformation (J3U)) - (|31|) to the initial data for the 
gravitating system of interest. In doing so, p becomes the proper rest-mass 
density and we substitute V^' = (c 2 /2)(l — /^) for the Newtonian potential 
uffl • (The coefficient is needed to ensure that ~ as A -> 0.) This 
defines a system (Sa) as the solution of this initial- value problem (Pa), and 
one expects that the solution fields admit expansions in powers of A, whose 
dominant terms have the same orders in A as in Eqs. (|3U|) - (j31|) . Thus, we 
introduce a conceptual family of gravitating systems, with the real system 
of interest being assumed to correspond to a given, small value Ao of the 
parameter A. It is then easy to check that, by adopting [M]a = A[M] and 
[T] A = [T]/V\ as the new units for the system (Sa) (where [M] and [T] are 
the starting units of mass and time), all fields become order A , and the 
small parameter A is proportional to 1/c 2 (indeed A = (c /c) 2 , where c is 
the velocity of light in the starting units). Hence, the derivation of the 1PN 
expansions and expanded equations is straightforward: all fields depend on 
the small parameter 1/c 2 , and one writes first-order expansions in this pa- 
rameter. However, there is still a crucial point to be precised, that will also 
play a role in the discussion of gravitational radiation. 
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3.1.2 The question of the time variable in the expansions 

The asymptotic expansions of the fields have first of all to be valid at fixed 
values of the time and space variables. Now, as we have just seen, the expan- 
sions are naturally formulated in units that depend on the small parameter 
A. Hence, it is not equivalent to assume expansions with x'° = T as the time 
variable, or with x° = cT, because the ratio x'°/x° — 1/c = A 1 / 2 /c depends 
on the small parameter. Note that, in the varying units, we have T = A 1/2 T 
where To is the true time (i.e. that measured in fixed units), hence cT = cqTq 
is (proportional to) the true time. But the true orbital velocities in the sys- 
tem Sa vary like A 1//2 , and the true orbital periods vary like A -1 / 2 . Hence it 
is T, not cT oc T , which remains nearly the same, as A is varied, for one 
orbital period of a given body in the Newtonian limit. Therefore, in this 
limit, thus for post-Newtonian expansions, one must take x'° = T as the time 
variable, not x° = cT. The latter choice corresponds to post-Minkowskian 
(PM) expansions, with wave equations for the gravitational potentials, which 
are relevant in the study of gravitational radiation (see Section 0J). 

3.1.3 Expansion of the metric and the field equation in the pre- 
ferred frame 

The leading expansion is that of the scalar field / = (3 2 (or rather, of the field 
V = (c 2 /2)(l — /), for this is really the field which is first-order expanded in 
1/c 2 below): 

/ = 1 - 2U/c 2 - 2A/c 4 + 0(l/c 6 ) (32) 
The space metric deduced from the Euclidean metric g° by Eq. (J3J) is 

9ij = g% + (2[//c>« + 0(l/c 4 ), h$ = (U^/ig^UM (33) 

(with g^- = g 0l i = 8^ in Cartesian coordinates), and the space-time metric 
immediately follows by P|): 

7oo = /, 7ij = ~9ij, loi = 0. (34) 

The mass-energy density a = T 00 may be written in the form 

a = a + ^i/c 2 + O(l/c 4 ), (35) 

where a = p is the conserved rest-mass density which is found at the first 
approximation. (Expanding, for a perfect fluid, the energy equation (J26|) . one 
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finds that p obeys the usual continuity equation; one also finds that mass 
is still conserved at the second (1PN) approximation.) The PN expansion of 
the field equation (jSJ) follows easily from Eqs. (J3*2"jl and (|33|) . Since the time 
variable is x'° = T for PN expansions, the time-derivative term in (jSJ) occurs 
with a c~ 2 factor as compared with the Af term, so that Poisson equations 
are obtained: 

AU = -AnGpo, (36) 
AA = -4ttG(t 1 + d 2 U/dT 2 . (37) 

3.2 Effects of a weak gravitational field on light rays 

3.2.1 Post-Newtonian equations of motion for a photon in the 
preferred frame 

The application to the situation of a light-like test particle in the field of a 
spherical body that undergoes a uniform translation in the frame E, with 
velocity dx/dT = V, is not difficult [6.. For a photon, the velocity u l is 0(c). 
Using Eqs. and (|3l?|). one gets the 1PN expansion of the coordinate 

acceleration (|18|) for a photon. In Cartesian coordinates for metric g°, this 
is 

% = U, - 2 {U„u>) § - ((<>),* + «>), - (<»)„) ^ + O . 

(38) 

The crucial point is that Eq. (|38jl . derived from Newton's second law, is 
nevertheless undistinguishable from the PN expansion of the equation for null 
space-time geodesies, as it is easily checked after Eq. (36)). Note that, 
in Eq. (|3B|) . the Newtonian gravity acceleration g = VqU (with components 
U t i in the Cartesian coordinates utilized) intervenes at the same order in A 
or 1/c 2 as the other terms (i.e. the order zero, although it is really a second- 
approximation formula: the first approximation gives simply du % jdT = 0). 

3.2.2 Transition to the moving frame 

Let us denote by E v the frame that undergoes a pure translation, with the 
assumed uniform velocity V of the massive body, with respect to the ether 
frame E. We may pass from E to Ev by a Lorentz transformation of the flat 
space-time metric 7 (that one whose line element is (ds ) 2 = c 2 dT 2 — dx l dx l 
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in Cartesian coordinates for the Euclidean space metric g°). At the 1PN 
approximation, a photon follows a null geodesic of the physical space-time 
metric 7, and, by the Lorentz transformation, the components of this space- 
time tensor transform like a (twice covariant) tensor, of course. This gives 
an easy way to get the PN equations of motion for a photon in the frame 
Ev- Neglecting 0(1/ c 4 ) terms, and except for an 0(1/ c 3 ) term in the 7^ 
components — such terms do not play any role in the 1PN expansion of the 
null geodesic equation |T5] — the components 7' after the Lorentz transform 
depend only on the "Newtonian" potential U, by just the same equations 
()32 )) -([34 j) as they do in the preferred frame. With the assumption of spherical 
symmetry, we have simply U = GM /r with M = J podV : then, Eqs. (132)) - 
(134)) are just the 1PN expansion of Schwarzschild's exterior metric. Whereas 
this spherical potential is, in general, not constant in the frame E, it is 
indeed constant in the frame Ev which moves with the spherical massive 
body (e.g. the Sun) that creates the relevant field, so that the geodesic 
equation is really that deduced from Schwarzschild's metric. We conclude 
that, to just the same level of approximation as in the 1PN approximation 
of GR, in particular neglecting the 7o« components of the metric, which are 
0(l/c 3 ), the present theory predicts exactly the same gravitational effects 
on light rays as the standard effects obtained in GR with the Schwarzschild 
metric — and this is true accounting for the preferred-frame effects, for they 
do not appear at this approximation. The details can be found in Ref. jU]. 

4 Gravitational radiation (isolated system) 

In GR and in the other relativistic theories of gravitation, most investiga- 
tions of gravitational radiation consider the so-called "weak-field approxima- 
tion", which consists of a linearization of the field equations. In particular, 
d'Alembert equations are obtained for the gravitational potentials, which is 
the basic reason why gravitational radiation then occurs in a tractable form. 
Obviously, this linearization can be done also in the present theory with its 
nonlinear wave-type equation (jSJ). 

4.1 Asymptotic post-Minkowskian approximation 

With the goal to develop mathematically clear approximation schemes, it 
seems desirable to find a definite asymptotic framework which would indeed 
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lead to a linearization of the field equations while keeping the propagating 
character. The PN scheme summarized in Subsection 13.11 does not fit with 
this aim, since it leads to the Poisson-type equations (jSfiJl and (|37jl . As 
it has been noted, this is due to the fact that, for consistency with the 
Newtonian limit, x'° = T must be taken as the time variable, in the varying 
units [M] A = A[M] and [T] A = [T]/VA, in which oc c 1 and all fields are 
order zero. It appears from ^3. 1.21 that, in order to have x° = cT as the time 
variable for different expansions, also based on this change of units, one must 
abandon the assumption that the velocity is low, and assume instead that it 
is order zero (in the starting units). 

4.1.1 Principle of the asymptotic PM approximation 

Thus, let us define a family (S' A ) of perfect-fluid gravitating systems by the 
following family of initial conditions (in fixed units): at the initial time, 

p ix) (x) = AV 1} (x) , (x) = Ap (1) (x) , (39) 

y (A) (x) = AV (1) (x), d T V w (x) = A9 r K (1) (x), (40) 

u (A) (x) = u (1) (x). (41) 

(Recall that V = (c 2 /2)(l — /).) This describes a limit in which the gravi- 
tational field becomes weak (due to a correspondingly diminishing density), 
but nevertheless the velocity u = |u| does not necessarily become negligible 
with respect to c — or, more exactly, u/c has some finite order of magnitude, 
which does not depend on A. We assume that, in a relevant time interval, 
the solution fields remain of the same order in A as in Eqs. (|39 p -()41| ) . Then, 
in the units [M] A = A[M] and [T] A = [T]/\/A, we have again vA oc c -1 , and 
all fields: p, p, V, dV/dx°, but the velocity u, are order zero: as to u, it is in 
fact order c. 

4.1.2 Zero-order PM approximation 

Since all fields but u are again order zero in the varying units, their (first- 
order) expansions have obviously the same form [Eqs. (JH2J)-(jHHJ)] as for the 
PN expansion. However, due to the fact that now u is 0(c), the zero-order 
term o"o in the expansion (|35jl of the active mass-energy density a = T 00 
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is not equal any more to the zero-order term p of the density of rest-mass 
in the preferred frame — as this is the case (201 for the PN approximation. 
Instead, <7o now contains kinetic energy in relativistic form: 

= Po (t«)o = Po (1 - Uq) = r, (42) 
where uo is the coefficient of c in the PM expansion of u: u = uqc + 0(c _1 ). 

Because x° = cT is now the time variable, the time-derivative term in 
Eq. (jSJ) for the scalar field / now does not contain the c~ 2 factor, hence it is 
expanded to a "flat space-time wave equation" with d'Alembert operator: 

DU = L/,0,0 - AU = AttGt (x° = cT). (43) 

The dynamical equations as applied to a perfect fluid, take the 

exact form 20J: 

9t (V/) + dj (iPfu j ) = (a/2) d T f + c~ 2 d T p, (44) 



d T + dj (imV) + T) k im j u h + t^u k = ijifg - g ij p :V (45) 

where 

+ ti = ^d T g jk . (46) 

When the PN approximation is used (time variable T, all fields of order zero), 
the zero-order expansion of these equations gives [20] the continuity equation 
and the classical fluid-dynamics equation in the presence of the gravity accel- 
eration £/j. (This equation can be put in the form of Euler's equation using 
the continuity equation.) In contrast, with the PM approximation (time 
variable x° = cT, and u of order c), we get a "continuity equation" which 
involves the density r (|4*2^) instead of the 0-order rest-mass density: 

d r + dj(Tui) = (x° = cT), (47) 

plus an equation meaning that the principal part of the acceleration is zero: 

d (Tui)+d j (rui4) = 0. (48) 

It should not surprise us that, when we do not impose the Newtonian limit, 
we do not get Newtonian equations indeed! The "asymptotic" scheme of 
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PM approximation that we have introduced in £ 14.1.11 (seemingly for the first 
time) is more general than the asymptotic PN scheme recalled in &3.1.11 in 
the sense that there is just more freedom left to the magnitude of the velocity. 
Thus, what happens here is that, due to the fact that u is order c, the acceler- 
ation U t i plays a role only at the following order (1PM approximation) — but 
the latter seems to be rather complicated. However, the zero-order PM ap- 
proximation is enough to assess the principal part of the radiative energy 
loss. 



4.2 0PM approximation of the radiative energy loss 

Considering an isolated self-gravitating system, assumed to be suitable for 
the PM scheme above (thus the gravitational field is weak, but the orbital 
velocities are not necessarily small), we wish to investigate the asymptotic 
behaviour, as Rq — > oo, of the energy variation (|2"8|) in the ball Q : |X| < Rq. 
More exactly, in the framework of the PM approximation, we investigate 
the limit, as Rq — > oo, of the principal part of dE/dT as A — > 0. Thus, we 
evaluate the flux $ (|2"Hj) at the zero-order PM (0PM) approximation. We 
have then from (J2J) and 

$ = $ (1 + 0(c- 2 )) , $ = -Zl /" d T U U k n k dS. (49) 

47r( ^ Jan 

4.2.1 Expansion of the retarded potential 

The relevant solution of Eq. (|4*H|) is the retarded potential: 

U(X, x°) = J Gr(X 'j g xI^~ X ' ) ^(x). (50) 

(The integral extends over the whole "retarded system", which is spatially 
bounded as is the system itself. Moreover, the point X shall be at a large 
distance from the retarded system.) The expansion of U as Rq = |X| — > oo is 
got by a standard procedure (see e.g. Stephani [21], pp. 125-126), which we 
try to make mathematically explicit — though without giving a convergence 
proof. One writes first the Taylor series expansion at x° — Rq of the retarded 
density: 

5R 2 

r(x, x°-R) = r(x, x° - R ) + 5R <9 r(x, x° - R ) + ^-«9 2 r(x, x°-Ro) + ..., 

(51) 
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where 

R = |X - x| , R = |X| , SR = R - R. (52) 
Inserting into (j51|) the first-order Taylor expansion of R at X: 

s «=% +o {-k)- (53) 

one gets: 

r(x,x°-i?) = t(x,x° - R ) + d T(x,x° - R ) 

+ »M-w (£y) + „.; (.) >(i) | 

where the points of suspension indicate omitted terms of order R® coming 
from the remainder in the infinite series (j51j) . Using then the zero-order 
expansion of 1/ |X — x| at X: 

1 l _ + (l 



IX-xl (55) 



together with (|H1|) . and integrating, gives the expansion of C/ up to l/Ro- 
U(X, T) M X l d\ X { Xi •■ „ ( 1 \ 

where overdot means partial derivative with respect to T, the points of sus- 
pension indicate omitted terms of order 1/Rq coming from time derivatives 
of order > 3 in (j51|) . and with 
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J rdV, di(X,T) = J zV(x,cT- \X\)dV(x), (57) 
JijiX, T) = J zVt(x, cT - |X|)dV(x). (58) 



4.2.2 Limiting value of the flux 

When taking the partial space-derivatives, we get from the definition (J57j) : 

ddi -X k 



dX k cR 
19 



d h (59) 



and the like for J^. We obtain therefore from (|56|): 



drUdkU 
G 2 



.V'.V y 

+ —r^J. 



cRl 



2c 2 R 3 13 



-X l X k d\ -X l X m X k 



C 2 ^ 



2c 3 ^ 



(We have dM/dT = by Eq. If37jl.) But it is also true that 



dt = 0. 



(60) 



(61) 



Indeed, the retardation being uniform in the definition ([57)1 . this is a stan- 
dard exact consequence of Eqs. (}4"T|) and (J4"%j) . for the continuum with density 
r = o" an d that moves with velocity cuo — the acceleration field of this con- 
tinuum being zero in view of Eq. (J48|) . On the other hand, in the next 
subsection we shall consider the same physical system at the Newtonian ap- 
proximation, thus equating r in the definition (JoTj) with the zero-order PN 
rest-mass density (which is consistent with ()42|) insofar as the actual velocity 
u is in fact negligible with respect to c). But at this approximation also it 
is true that di = 0, because the total momentum of an isolated Newtonian 
system is conserved. 



Thus, there is no dipole term in fact. From (J4T?|) . ()60|) and (j61|) . we get 

-G 2 /•/••■ • a 2 



-AnG lim $ = r- 

i?0^oo u 4c 5 



(62) 



where duo is the solid angle element and n l = X 1 /Rq. (Note that tensor J is 
constant on the sphere |X| = Rq.) Using the integration formula for n l n : 'n k n l 
21 (Eq. (90.20)), it follows that 



lim 



G 
60c5 



2tr J + tr J 



(63) 



4.3 Peters-Mathews coefficients for a binary system 

The energy loss E = — ($o)iim is the principal part of the variation of the 
total energy of our isolated weakly-gravitating system, due to its emission 
of gravitational waves. As it was already mentioned at several places in 
this Section, this unique system may consistently be treated in both the PM 
approximation and the PN one. The (more general) PM approximation has 
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allowed us to calculate an approximate value of the energy loss, and now we 
are going to use the OPN (Newtonian) approximation to see how the energy 
loss changes the relative motion in the first approximation. To do this means 
assuming: i) that the OPN approximation of the energy of the system (the 
Newtonian energy, which is an exact constant in the framework of the OPN 
approximation), actually varies as the exact energy of the system — and ii) 
that the latter variation is well approximated by the 0PM calculation. 6 We 
consider more particularly a binary system, made of two massive bodies that 
are far enough from each other, so that tidal interactions are negligible and 
each body has a rigid motion. We assume, moreover, that the self-rotation 
of each body is not modified by the energy loss, which thus affects merely 
the relative motion of the two bodies, whose motion may then be analysed 
as a Newtonian 2-body problem for two point masses mi and m 2 , located at 
xi(T) and X2(T). However, since E has been evaluated in the preferred frame 
E, this Newtonian problem is formulated in that frame (which is indeed an 
inertial frame for the OPN approximation of the theory). Setting 

x = xi - x 2 , u = x, (64) 

we have thus 

ii = — fcx/r 3 , k = G{m\ + m 2 ), r = |x| . (65) 

The quadrupole tensor is then: 

J = miXi <g> xi + m 2 x 2 <g) x 2 . (66) 

Introducing the mass center a(T), such that miXx + m 2 x 2 = Ma, where 
M = mi + m 2 , we obtain 

J = J + Ma ® a, J = m 1 (x 1 — a) <g> (x x — a) +m 2 (x 2 — a) <g> (x 2 — a). (67) 

Since the system is isolated, the velocity V = a is a constant. It follows that 

J = J. (68) 

6 We note that the same two assumptions have to be done when one applies the 
quadrupole formula of GR, and that assumption i) may seem strong. A priori, one might 
expect that the energy loss found at the 0PM approximation would affect also the higher- 
order terms (possibly even only the higher-order terms!) in the PN energy. Of course, the 
Newtonian energy dominates over the PN corrections to it. Assumption i) is valid if this 
remains true for the corresponding rates, which should be true if the PN expansions are 
uniform with respect to time. 
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We get from 

J = /xx ® x, /x = m,\m2jM. 
Three time differentiations of gi ye us with and 



[3(u.n) n <g> n — 2(n (g> u + u (g> n)] , n = — 



whence follows: 
tr J 



2kfj. 



4k 2 fi 2 



-u.n, 



tr y 



[8u 2 -7 (u.n) 2 ] . 



Using Eq. (|63|). this leads finally to 

8 G k 2 fi 2 



OJhm 



15 c 5 



2u 2 - 



u.n 



i.e. the so-called Peters-Mathews coefficients (see e.g. Will [IB], Eq. 



^ 5 ^dipole 



0. 



(69) 
(70) 
(71) 

(72) 
10.84)): 

(73) 



There is no average to be performed over one binary's orbit in Eq. (p)3*|) nor in 
Eq. (|72j) — in contrast with what is often found in the literature. If we start 
from the "quadrupole formula of GR" (the equivalent, for GR in harmonic 
coordinates and for the RTG, of Eq. (jfiH)) here: cf. Ref. [13], Eq. (110,16), 
and Ref. [TTj, Eq. (15.71), where there is indeed no time average): 



E 



G 



45c 5 



tr D 



itr 2 J 

3 



D 



- (trJ)I 

3 V ; 



(74) 



[12 u 2 - 11 (u.n) 2 ] 



(75) 



and apply Eq. (fTTj) . we find 

^ _ 8 G k 2 fi 2 

which is indeed the Peters-Mathews formula of GR in harmonic coordinates 
[TT)] and of the RTG jllj . without any time average (in contrast with Ref. [TH] . 
Eq. (10.80)). The time-average may appear physically justified, for one ex- 
pects that the time-period of the gravitational waves should be of the order of 
the orbital period. But, as we have outlined at the beginning of this subsec- 
tion, it seems in fact reasonable to admit that the rate of the radiative energy 
loss equals the decrease-rate of the Newtonian energy. Once this equality is 
admitted, it implies a "Peters-Mathews formula" without any time-average. 



22 



5 Conclusion 



We have given a nearly self-consistent summary of the "scalar ether-theory" : 
Eqs. d2J for the gravity acceleration, (J3J) and (J3J) for the metric, and (0) for 
the scalar field, can be taken axiomatically. We have more detailed its dy- 
namics (Eq. (j!5)) for a test particle, Eq. (|2U|) for a continuum). It should be 
noted that this dynamics is independent of the assumptions made as to the 
form of the metric. A detailed motivation for the theory, together with a 
summary of its construction, can be found in Refs. [T2| or [H|. 

This scalar theory gives the same predictions for light rays as the stan- 
dard predictions deduced from Schwarzschild's metric in GR. This is be- 
cause the theory predicts Schwarzschild's metric (Eq. ©), and because it 
turns out that the preferred-frame effects do not play any role at the first 
post-Newtonian (1PN) approximation for photons. In turn, the fact that 
this central metric of GR also does occur in this scalar theory, is related to 
its preferred-frame character, which allows to postulate gravitational time- 
dilation and space-contraction in a way that is not space-time covariant. 

To investigate gravitational radiation, we have introduced an "asymp- 
totic" form of the post-Minkowskian (PM) approximation, based on an ex- 
plicit family of initial-value problems for a perfect-fluid system (Eqs. (I39j) - 
(141(0 . This seems to be new. We did not prove the existence of the solutions, 
nor did we prove the truly asymptotic character of the expansions, but we 
consider it likely that both could be proved (in suitable function spaces). 
The important mathematical work [221, that is devoted to GR (whose gravi- 
tational field equations are much more complex), does give such justifications, 
but it is restricted to the vacuum Einstein equations. As it was the case for 
our "asymptotic" PN scheme [201 > the present PM scheme justifies, by invok- 
ing a change of units, to take 1/c 2 as the small parameter, while considering 
that (nearly) all fields are order zero — but all fields depend on the small 
parameter and hence must be expanded. The difference between these PN 
and PM schemes is that, in the latter, the time variable is x° = cT, not T, 
and the velocity is 0(c) only. But the latter is true in the changed units 
and in an asymptotic sense (i.e. for the conceptual one-parameter family 
of gravitational systems), physically it means that the velocity field can be 
non-negligible with respect to c, but it does not have to be so. Thus this 
PM scheme is just more general than the corresponding PN scheme. This 
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is the reason why the energy rate calculated in the 0PM approximation can 
be used to investigate the time evolution of the system, when the latter is 
considered at the PN approximation. Thus, within this framework, it seems 
that one can reasonably assume that it is not undispensable, in a first ap- 
proximation of the "radiation reaction effects" , to make a complete analysis 
of the relativistic 2-body problem up to the order 1/c 5 or higher. Such a 
complete analysis has been investigated for GR (cf. Damour j21] and refer- 
ences therein), but it is very complex. 

Using the PM approximation, we have shown that, in this scalar theory 
with a preferred reference frame, there is a flux of gravitational energy which 
emanates from an isolated weakly-gravitating system, and which goes to in- 
finity; that this limit flux is positive, Eq. (J53j) (hence the system loses energy); 
and that it starts only from quadrupole terms, the monopole and dipole terms 
being zero. When equated to the energy rate of a binary system in the New- 
tonian approximation, this energy loss gives a Peters-Mathews formula with 
well-ordered positive coefficients, again without any dipole (Eq. (|72)l). Since 
these coefficients differ from those in GR or in the RTG, a different value of 
the energy loss would be obtained (albeit of the same order of magnitude), 
if one would apply Eq. (j72j) while otherwise keeping the parameters of the 
2-body system at the values found [22] from fitting the pulse data with a 
"timing model" based on GR. However, there is essentially only one input 
data for the timing model, namely the list of the time intervals between the 
successive pulses emitted by the pulsar (in short "the pulse data"), and if the 
theory is changed, then, after fitting the pulse data by a timing model based 
on the new theory, different values will be found also for the parameters of 
the 2-body system, e.g. for the eccentricity and for the masses of the two 
bodies. Therefore, the author feels it plausible that the pulse data of binary 
pulsars might be nicely fitted by building a timing model entirely based on 
this theory, in the same way as they have been fitted in GR [2*3] . What can 
only be stated for now, is that the theory predicts an energy loss with the 
same order of magnitude to that which is found in GR and which fits with 
the observations. As emphasized by Will jTH], most alternative theories do 
not predict the same order of magnitude as in GR for the energy rate, and 
not rarely they even predict the wrong sign — e.g. this is the case for Rosen's 
bimetric theory that looks much closer to GR than looks the present 
scalar theory. 
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